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ABSTRACT 

We show that the relativistic two-stream instability can naturally generate strong magnetic fields with 
10~^-10^^ of the equipartition energy density, in the collisionless shocks of Gamma- Ray-Burst (GRB) 
sources. The generated fields are parallel to the shock front and fluctuate on the very short scale of the 
plasma skin depth. The synchrotron radiation emitted from the limb-brightened source image is linearly 
polarized in the radial direction relative to the source center. Although the net polarization vanishes 
under circular symmetry, GRB sources should exhibit polarization scintillations as their radio afterglow 
radiation gets scattered by the Galactic interstellar medium. Detection of polarization scintillations 
could therefore test the above mechanism for magnetic field generation. 

Subject headings: gamma rays: bursts — magnetic fields — instabilities 



1. INTRODUCTION 



Narayan, & Piran 1996) 



Cosmological 7-ray bursts (GRB) are believed to be pro- 
duced in the fireballs of very energetic explosions, when a 
large amount of energy, E ~ gj.g^ jg released over a 

few seconds in a small volume in with a negligible baryonic 
load, Mc^ <C E (see Piran 1999 for a review). Most of the 
energy is eventually transferred to the baryons which are 
accelerated to ultra-relativistic velo cities with a Lorentz 
factor 7 ^ E/M c^ - lO^-lO^ (e.g., |Shemi fc Piran 1990; 
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eth ' 



eth 



(1) 



Here Ub = B^/Stt and Ue are the magnetic and electron 



energy densities and eth 



pC (7p — 1) is the total ther- 



Paczynski 199Cl| ). A substantial fraction of the kinetic en- 



ergy of the baryons is transferred to a non-thermal popula- 
tion of rela tivistic electrons throug h Fermi acceleration at 
the shock (Meszaros & Recs 1993). The accelerated elec- 



trons cool via inverse Compton scattering and synchrotron 
emission in the post-shock magnetic fields and produce 



mal energy density behind the shock; where nip is the pro- 
ton mass, n is the proton number density, and 7^ is the 
mean thermal Lorentz factor of the protons. The observed 
afterglow spectra and lightcurves typically yield values of 
the magnetic energy parameter ranging from cb ~ 0.1 
(Waxman 1997; Wijers & Galama 1998), down to 10"^ 
(Gra not et al. 1998) or ev en es ~ 10^^ (Galama et al. 



Katz 1994; 


Sari, Narayan, & Piran 1996; Vietri 1997 




Wax- 


man 1997a 


; Wijers, Rees, & Meszaros 1997). The shock 



could be either internal due to collisions between fireball 
shells caused by source variability (Paczyhski & Xu 1994; 



Rees fc Meszaros 1994 ) , or external due to the interaction 
of the fireball with the surro unding interstellar medium 
(ISM; ^ eszaros fc Rees 19931 ). The radiation from inter- 
nal shocks can explain the spectra (Pilla & Loeb 1998) 
and the fast irregular variability of GRBs (Sari & Piran 
1997a), while the synchrotron emission from the external 
i shocks provides a successful model for the broken powcr- 



1999; Vreeswijk et al. 1999) — all below the equipartition 
limit ~ 1. 

The existence of strong magnetic fields is naturally ex- 
pected in the compact environments of potential GRB pro- 
genitors. First, the field might originate from a highly 
magnetized stellar remnant, such as a neutron star, with 
B < 10^^ gauss. Second, a turbulent magnetic dynamo 
could amplify a relatively weak seed magnetic field in the 
vicinity of the progenitor. This process, however, requires 
the turbulence to be anisotropic and have a nonzero to- 
tal helicity, v • (V x v) 7^ 0. A similar mechanism, called 
the Q -r^-dynamo, mig ht opera te in rapidly rotating ob- 



Inw opootpn and omooth tomporol bohiwioii of aftorgloTi 



jccts (Thompson 1994, see also Meszaros, Laguna, fc Rees 
1993[). Finally, the magnetic shearing instability ( Palbu! 



(e.g., Waxman 1997a, b). In both caeoe, strong magnetic 



& Hawley 1991) could amplify the magnetic field (but not 
the flux) in strongly sheared flows. The e- folding time 
for this instability is approximately the rotation period, 
which decreases with radius as due to angular mo- 
mentum conservation in the outflowing wind. Thus, being 
possi bly important in the early stages of t he fireball expan- 
sion (Narayan, Paczyhski, & Piran 1992), this instability 
is inefficient at large radii, where its e-folding time greatly 
exceeds the dynamical time of the fireball. 

In contrast to such progenitor environments where large 
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fields are required behind the shocks at all times in order 
to fit the observational data. 

The properties of the synchrotron emission from GRB 
shocks are determined by the magnetic field strength, B, 
and the electron energy distribution behind the shock. 
Both of these quantities are difficult to estimate from first 
principles, and so the following dimensionless parameters 
are often used to incorporate modeling uncertainties ( Sari J 
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magnetic fields are natural, there is currently no satis- 
factory explanation for the origin of the strong magnetic 
fields required in GRB shocks (see discuss i ons in Thomp-| 
son 1 $94|; Meszaros, La guna, fc Rees 1993| ; pari, NarayanJ 



Meszaros, Laguna, fc Rees 1993 ; 3ari, Narayan 
1996| ). Compression of the ISM magnetic field 
in external shocks yields a field amplitu de B ~ 7-Bt.sm 



10 "^( 7/10^ ) gauss, which is too weak ( sari, Narayan, fc 
Piran 1996| ) compared to the required equipartition value 



Beq ^ 50(7/10^)(nigM/l cm~'^)^/^ gauss, and can ac- 
count only for eb — {B/B^qY ^ 10"^^. Alternatively, 
some magnetic flux might originate at the GRB progen- 
itor and be carried by the outflowing fireball plasma (or 
by a precursor wind). Because of flux freezing, the field 
amplitude would decrease as the wind expands. In this 
case, only a progenitor with a rather strong magnetic field 
^ 10^^ gauss might produce sufficiently strong fields dur- 
ing the GRB emission. However, since the field amplitude 
scales as B (X V~'^^^ for an expanding shell of volume V, 
even a highly magnetized plasma at i? '--^ 10^ cm would 
possess only a negligible field amplitude of ^ lO^'^ gauss, 
or cb <, 10"'', at a radius of R ^ 10^^ cm, where the af- 



terglo w radia tion is emitted'^ (see also Meszaros, Laguna 



fc Rees 1993). Moreover, the emitting material behind the 
external shock is continuously replenished by the ISM, and 
so the field originally carried by the fireball ejecta cannot 
account for the afterglow radiation. 

None of the above mechanisms is capable of generating 
near equipartition magnetic fields in the external shocks 
which produce the delayed afterglow emission. In this pa- 
per, we propose a different, universal mechanism of mag- 
netic field generation in GRB shocks. It involves the rela- 
tivistic generalization of the two-stream (Weibel) instabil- 
ity in a plasma. This instability is driven by the anisotropy 
of the Particle Distribution Function (PDF) and, hence, 
could operate in both internal and external shocks. Our 
main results are as follows: 

1. The characteristic e-folding time in the shock frame 
for the instability is ^ 10"'' s for internal shocks and 
10^"* s for external shocks. This time is much shorter 
than the dynamical time of GRB fireballs. 

2. The generated magnetic field is randomly oriented 
in space, but always lies in the plane of the shock 
front. 

3. The instability is powerful. It saturates only by non- 
linear effects when the magnetic field amplitude ap- 
proaches equipartition with the electrons (and pos- 
sibly with the ions). 

4. The instability isotropizes the PDF and, thus, effec- 
tively heats the electrons and protons. 

5. The characteristic coherence scale of the generated 
magnetic field is of the order of the relativistic 
skin depth, i.e. ~ 10'^ cm for internal shocks and 
~ 10^ cm for external shocks. This scale is much 
smaller than the spatial scale of the source. 

6. The mean free path for Coulomb collisions is larger 
than the fireball size. However, the randomness of 

^Both the magnetic field energy density and the thermal energy 
shocks are generated, the plasma is heated due to the dissipation of 
far below equipartition in the post shock gas. 



the generated magnetic field provides effective colli- 
sions due to pitch angle scattering of the particles in 
an otherwise collisionless plasma and, thus, justifies 
the use of the magncto-hydrodynamic (MHD) ap- 
proximation for GRB shocks. The magnetic fields 
communicate the momentum and pressure of the 
outflowing fireball plasma to the ambient medium 
and define the shock boundary. 

The above mechanism results in tangential magnetic 
fields near the apparent limb of the source. Hence, the 
long-term synchrotron emission from the limb would be 
linearly polarized along the radial direction relative to the 
source center. Although the net polarization of a circularly 
symmetric source is zero, scattering of the radio afterglow 
emission of GRBs by the intervening Galactic interstellar 
medium would break the symmetry in the source image 
and result in polarization scintillations. This effect can be 
used to test the reality of our proposed mechanism for the 
generation of magnetic fields in GRB blast waves. 

The outline of the paper is as follows. The physical 
mechanism of the instability is discussed in The gen- 
eration of magnetic fields in internal and external shocks 
is discussed in §1]. In §^ we predict the polarization scin- 
tillation signal in our model. Finally, §g| summarizes our 
main conclusions. 

2. THE TWO-STREAM INSTABILITY 

The instability under consideration was first predicted 
by Weibel (1959) for a non-relativistic plasma with an 
anisotropic distribution function. The simple physical in- 
terpretation provided later by Fried (1959) treated the 
PDF anisotropy more generally as a two-stream configu- 
ration of a cold plasma. Below we give a brief, qualitative 
description of this two-stream magnetic instability. 

Let us consider, for simplicity, the dynamics of the elec- 
trons only, and assume that the protons are at rest and 
provide global charge neutrality. The electrons are as- 
sumed to move along the a;-axis (as illustrated in Figure 
|l]) with a velocity v = zLxv^ and equal particle fluxes in op- 
posite directions along the x-axis (so that the net current 
is zero). Next, we add an infinitesimal magnetic field fluc- 
tuation, B = zBz cos{ky). The Lorentz force, — x B, 
deflects the electron trajectories as shown by the dashed 
lines in Figure [|. As a result, the electrons moving to 
the right will concentrate in layer I, and those moving to 
the left - in layer II. Thus, current sheaths form which 
appear to increase the initial magnetic field fluctuation. 
The growth rate is F = UpVy/c, where cijp — (47re^ n/m) is 
the non-relativistic plasma frequency (Fried 1959). Simi- 
lar considerations imply that perpendicular electron mo- 
tions along y-axis, result in oppositely directed currents 
which suppress the instability. The particle motions along 
z are insignificant as they are unaffected by the magnetic 
field. Thus, the instability is indeed driven by the PDF 
anisotropy and should quench for the isotropic case. 

The Lorentz force deflection of particle orbits increases 
as the magnetic field perturbation grows in amplitude. 
The amplified magnetic field is random in the plane per- 
pendicular to the particle motion, since it is generated 

of the fireball scale as oc V"*/-^ for adiabatic expansion. However, when 
the fireball kinetic energy, and the magnetic energy parameter decreases 
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from a random seed field. Thus, the Lorentz deflections 
result in a pitch angle scattering which makes the PDF 
isotropic. If one starts from a strong anisotropy, so that 
the thermal spread is much smaller than the particle bulk 
velocity, the particles will eventually isotropize and the 
thermal energy associated with their random motions will 
be equal to their initial directed kinetic energy. This final 
state will bring the instability to saturation. 

We note the following points about the nature of the 
instability: 

1. The instability is aperiodic, i.e., Rew = 0. Thus, it 
can be saturated by nonlinear effects only, and not 
by kinetic effects such as collisionless damping or res- 
onance broadening. Hence, the magnetic field can be 
amplified up to high values. 

2. Despite its intrinsically kinetic nature, the instabil- 
ity is non- resonant,* i.e., it is impossible to single out 
a group of particles that is responsible for the insta- 
bility. Since the bulk of the plasma participates in 
the process, the energy transferred to the magnetic 
field could be comparable to the total kinetic energy 
of the plasma. Hence, the instability is powerful. 

3. The instability is self-saturating. It continues un- 
til all the free energy due to the PDF anisotropy is 
transferred to the magnetic field energy. 

4. The generated magnetic field always lies in the plane 
perpendicular to the initial anisotropy axis of the 
PDF, i.e., to the shock propagation direction. 

5. The produced magnetic field is randomly oriented 
in the shock plane. The Lorentz forces randomizes 
particle motion over the pitch angle and, hence, in- 
troduces an effective scattering process into the oth- 
erwise collisionless system. This validates the use of 
the MHD approximation in the study of collisionless 
GRB shocks. 



Sagdeev fc Galeev (1969) and [Moiseev fc Sagdeev (1963)| 
provide a kinetic, non-relativistic treatment of the insta- 
bility in both the linear and the quasi-linear regimes, and 
apply the theory of collisionless shocks to space plasmas. 
In the next section we will extend their analysis to the case 
of ultra-relativistic GRB shocks. 

3. MAGNETIC FIELD GENERATION IN GRB SHOCKS 

We consider a GRB shock front expanding at a Lorentz 
factor, 7sh, behind which the particles have a thermal 
Lorentz factor, 7. In this section, we will derive equa- 
tions which are equally applicable to electrons and pro- 
tons, whichever species dominates the growth of the insta- 
bility. Later, we shall use the subscripts "e" and "p" to 
denote electrons and protons, respectively. We calculate 
all quantities in the comoving frame of the shock. 

A fully kinetic, relativistic treatment of the magnetic 
two-stream (Weibel) instability is a complicated task. The 
dispersion relation for a simplified "water-bag" PDF was 
derived by Yoon & Davidson (1987) and is given by equa- 
tion ( A_2) in Appendix This dispersion relation implies 
that only a range of modes above a critical wavelength will 



grow [cf. Eq. (A3)]. Naturally, the mode with the largest 
growth rate, Fmax, dominates and sets the characteristic 
length-scale of the magnetic field fluctuations, A ~ ^mix- 
The ultra-relativisti c ex pressions for Fmax and fcmax are 
given by equation (A5) for a strong initial anisotropy. 
We write the corresponding e-folding time and correlation 
length of the fleld as 



1/2 

r ~ , A 



1/2 



(2) 



We can now estimate the nonlinear saturation ampli- 
tude of the magnetic fleld. The instability is due to the free 
streaming of particles. As the field amplitude grows, the 
transverse defiection of particles gets stronger, and their 
free streaming across the field lines is suppressed. The typ- 
ical curvature scale for the deflections is the Larmor radius, 
P = v±b/^c — (lis ~ ly^'^mc^/eB, where v±b and 7j__b 
are the transverse velocity and Lorentz factor of a particle 
relative to the local magnetic fleld. On scales larger than p, 
particles can only move along fleld lines. Hence, when the 
growing magnetic flelds become such that /cmaxP ~ 1, the 
particles are magnetically trapped and can no longer am- 
plify the field. Assuming an isotropic particle distribution 
at saturation (7_l_b ~ 7), this condition can be re- written 
as 



(7 + 1) 

7T!,c^n(7 — 1) 2\/27 



(3) 



For 7^1, this corresponds to a magnetic energy density 
close to equipartion with the amplifying particles. Inter- 
estingly, one may obtain the same result following a dif- 
ferent analysis. First, the instability leads to a growth 
of the field amplitude [as given by the last term in Eq. 
(Al), ~ v-9x/]- Second, nonlinearity leads to the trans- 
fer of energy to shorter wave lengths, k > fccHt, where 
the fluctuations are damped [as described by the second 
term in Eq. (Al), ^ (e/c)v x B-9p/]. Thus, the steady 
value of B is determined by balancing these two pro- 
cesses. Equating these two terms and replacing dx by 



^crit — ^max /V2^Pp/V2, yields 

v±Bkci-itf ~ v^BeBf/mc^j^B- 



(4) 



The field strength estimated here is equivalent to that 
given in equation (^ to within a factor of order unity. 

Direct computer simulations of the instability in both 
non-relativistic and relativistic electron plasmas confirm 
that the saturation occurs at slightly sub- equipartition val- 
ues of i? (see, e.g.. palifano ct al. 19*98: Kazimura et al. 



1998; 



yang et al. 1994| ; |Wallace fc Epperlein 199l|) 



BySn 



mc^nl^ — 1) 



ry- 0.01 -0.1. 



(5) 



where we introduced the efficiency factor 77 <^ 0.1. The pre- 
cise saturation level depends on the nonlinear modification 
of the PDF during the instability which is not accounted 
for by our linear analysis. We shall retain the efficiency 
factor, f], in our estimates. 



This instability may be treated as an analog of the fire-hose instability in the absence of the external magnetic field. 
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Note that the thermal Lorentz factor of particles, 7, 
varies in time as the instability develops. Due to par- 
ticle scattering by the generated magnetic fields, an ini- 
tially highly anisotropic PDF with 7 <C 7sh will eventu- 
ally evolve to an isotropic, ring-like distribution, for which 
7 ~ 7sh- Thus, the spatial scale and amplitude of the re- 
sultant magnetic field, given by equations (||) and (||), will 
evolve during the lifetime of the instability because they 
are functions of 7. In estimating these values at a GRB 
shock when the instability saturates, we take 7 ~ 7sh- The 
e-folding time for the instability is independent of 7 in the 
case of strong anisotropy. In the case of weak anisotropy, 
7 ~ 7sh, the "water-bag" model used here is formally 
invalid, but the comparison of the ultra-relativistic re- 
sults with non-relativistic results (e.g., Moiseev & Sagdeev 
1963) suggests that the instability quenches and the e- 
folding time scales as 



A/ccx [(e|| -e_L)/e||] 



-3/2 



(6) 



where ey and e± are the average energies of particle mo- 
tions along the direction of shock propagation and trans- 
verse to it. The field correlation length follows a similar 
scaling. 

The diffusive decay time of the generated magnetic field 
is Tdiff ~ l/rys/c^axi where rjB = mc^ i/coii/47rne^ is the 
magnetic diffusivity and ;/coii is the particle collision fre- 
quency. Hence, the diffusion time 



Tdiff 



7/ 



t'coU, 



(7) 



is much longer than the fireball expansion time since the 
particle collision frequency in the fireball plasma is negli- 
gible. Thus, the magnetic field is not expected to dissipate 



its ent!igy Oluuically uvei Lhe fiiuball lifuLime. NuLu LhaL 
magn etic fields cannot be produced during the optically - 
thick phase of the fireball, because Compton scattering on 
the photons rapidly removes any anisotropy of the PDF. 
Next, we consider two types of GRB shocks in which mag- 
netic fields might be generated. 

3.1. Internal Shocks due to Shell Collisions Inside the 

Fireball 



Rajbid variability of a GRB source results in a fireball 



ojp^ = 9.0 X lO^ni/2 and Wp^ = 2.1 x lOV/^ s'^ 

where n is in cm^'^. 

For simplicity, we consider the collision of two identical 
shells. In the center of mass frame, the interaction of these 
collisionless shells yields a state of two inter-penetrating 
plasma streams, which is readily unstable to the genera- 
tion of magnetic fields. Since 3> ^pp, the instability 
grows faster for the electrons than for the protons and 
so the electrons dominate the magnetic field generation 
process at early times. The electron instability saturates 
when the magnetic energy density becomes comparable 
to the electron energy density, ymtnmeC^ . This energy is 
still much smaller than that associated with the protons^. 
Thus, when the instability saturates for the electrons, it 
could still continue on a longer time-scale for the protons. 
The protons dominate energetically and could lead to near 
equipartition magnetic energy with 



es^r^-^O.l. 



(8) 



From equation (g) we get the characteristic scale length 
and growth time of the instability for the protons. 



A~ 2 X 10^ 



r ~ 6 X 10" 



-1/2 



/7inty/2 

3 X 1010 cni-3 ; V 4 / 

-1/2 

/7int\ 

3 X IQio cm-3 J V 4 / 



cm, (9a) 



1/2 



s. (9b) 



These quantities are decreased by a factor of 
{nip/raeY/'^ — 43 for the electrons. 

The generation of magnetic fields in counter-streaming, 
electron-positron plasmas has been extensivel y studied nu- 



merically using particle-in-cell codes 

1998| ; ICalifano et al. 199S| ; [Ymg et al 



al 



Ka zimur a et 
1994|) 



which IS composed ot thm layers (^shells) moving with 
different Lorentz factors. To produce the observed non- 
thermal 7-ray spectrum, the shells must collide at suf- 
ficiently large radii where the internal shock region is 
optically-thin to both Compton scattering and e^e^-pair 
production. The collision should also occur before the fire- 
ball slows-down on the ambient medium. These conditions 
imply that the internal shock be mildly relativistic, with 
a Lorentz factor 7int of order a few in the center of mass 
frame of the colliding shells (see Piran 1999 for more de- 
tails). Prior to a collision, the electrons and protons in 
the colliding shells are cold relative to their bulk Lorentz 
factor, 7e.p < 7int- As typical parameters for the shells we 
assume a plasma density of n w 3 x lO^" cm^'"^, 7int — 4, 
and initial thermal Lorentz fa ctors 7p^e ~ 2 (see e.g.. Pi- 
ran 1999; Pilla fc Loeb 199^ ). The plasma frequencies 
for the electrons and protons are given by the relations, 

^We assume that the dominant ion species in the relativistic GRB 
species is straightforward. 



A 

clear visual demonstration of the magnetic field amplifi- 
cation process is provided by Figure 2 of Kazimura et al. 
(1998). The rapid generation of a strong, small-scale mag- 
netic field occurs at the interface of the colliding streams, 
and is followed by the gradual modification of the field 
structure around the interface, due to the nonlinear sat- 
uration and relaxation of the particle velocity anisotropy. 
The inferred value of 77 ^ 0.01 — 0.1 is generic (Kazimura 
et al. 1998| Yang et al. 1994). The amplication process 



produces also random electric fields with an energy den- 
sity that is at most comparab le to that of the mag netic 
component, {E"^) ~ {v/c){B'^) (Kazimura et al. 1998). 

Unfortunately, no simulations were performed so far for 
colliding electron-proton plasmas. Numerical simulations 
of a plasma with species of somewhat different masses sug- 
gests that the energetics of the proce ss is indeed dominated 
by the heavier species (Arons 1996). Nevertheless, direct 
relativistic simulations with dynamical protons and elec- 
trons are required in order to assess the saturation am- 
plitude of the magnetic field in GRB shocks of different 
properties. 

If the colliding shells do not possess similar densities, 
then the growth rate of the instability decreases or even 
shuts off beyond a particular density contrast, as discussed 
in Appendix In this regime, the shock may be domi- 
nated by electrostatic (Langmuir) turbulence. Unless the 

wind is protons. The generalization of our discussion to heavier ion 



5 



outflowing plasma is already contaminated by strong mag- 
netic fields, the synchrotron emission from the collision 
of shells with very different densities would therefore be 
weak. 

3.2. External Shock due to the Interaction of the Fireball 

with the ISM 

Eventually, the fireball slows down due to its interaction 
with the surrounding ISM. The external shock produced 
by this interaction yields the delayed afterglow emission 
and is assumed to carry a strong magnetic field. As the 
shock propagates into the ISM, the fresh electrons and pro- 
tons are reflected from the magnetized shock front back 
into the ISM. Thus, a two-stream state forms in the co- 
moving frame of the shock and a magnetic fleld is amplifled 
in the ISM, just in front of the shock. 

We assume that the fraction of reflected particles is of 
order unity, and so the above two-stream state is analogous 
to that produced in internal shocks^. The instability flrst 
acts on the electrons. The correlation scale and saturation 
amplitude of the field are given by equations and (||). 
Since the magnetic field is generated upstream and then 
transported downstream, we need to take account of the 
compression factor at the shock. Given the jump condi- 
tions for a relativistic shock, we have A = A'/47sh (where 
"prime" denotes the parameter in front of the shock) , while 
the ratio of the magnetic to thermal energy remains con- 



Foi' ail ISM deiislLy Uism ~ 1 "^lu , we thei'cfoi'c 



stant. 
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<^Be^V [ ~ 5.5 X 10 



Ae = 3 X 10 



4 X 10" 



( Ish 

\ 10 

-4 fish 

10 



-1/2 



1/2 



1 cm^-^ 



-1/2 



cm. 



VI c 



ISM 



-1/2 



(10a) 

(10b) 
(10c) 



where rj i = (yy/O.l) and the subscript "e" denotes am- 
plification of the magnetic field by the electrons only. The 
magnetic energy parameter is still normalized relative to 
the proton thermal energy. 

When the instability of the electrons saturates, further 
amplification by the protons may become important. The 
magnetic field is amplified in a thin layer in front of the 
shock, the width of which is of order the Larmor ra- 
dius of the protons^ at the shock. The time available 
for field amplification is, thus, roughly the crossing time. 



amp 



3{mp/me)^^^{Bsi^t,e/B)Tp, where Bgat.e denotes the field 
stre ngth after saturation on the electrons [as given by 
Eq. (10a)]. On the other hand, the growth of the field to a 
sub-equipartition amplitude with the protons would take 

at least tgrowth = Tp ln(Bsat,p/Ssat,e) ~ Tp\n{mp/me)^/'^ ~ 

3.8Tp, i.e. comparable to the time available for the ampli- 
fication of i?sat,e up to Bgat.p = (™p/?7T-e) ^^^-Bsat,e ■ How- 

ever, since the growth of the field near saturation is slower 
than that during the linear stage of the process, the Weibel 
instability may not be able to build the field up to equipar- 
tition with the protons and yield es ~ 77. Whether max- 
imal amplification of the magnetic field can occur in this 
environment is uncertain'^ and can be found only through 
detailed numerical simulations. We thus conclude that 
the most robust prediction for the value of the magnetic 
field energy is ~ ri{me/mp), but somewhat higher val- 
ues are also possible, so that 



5 X 10 ^77.1 < es ^ O.I77.1. 



(11) 



The predicted range of eb ^ 10^^-10"^ matches the 
results from modeling of recent afterglow data. Wijers 
& Galama (1998) show that X-ray to radio spectrum 
of GRB970508 afterglow is consistent with the values of 
€3 ~ 0.07, indicating the proton dominated regime of field 
generation. The field energy density for GI IB990123 and 
GRB 971214 is estimated to be £3 ^ 10~^ ([Galama et af 
19991) , and for GRP980703 6 x 10"'"^ ([Vreeswijk eFal 



1999), consistently with the electron dominated regime 



4. POLARIZATION SCINTILLATIONS 

In the previous section, we have found that the char- 
acteristic time it takes the magnetic field to grow up to 
equipartition values is orders of magnitude shorter than 
the dynamical time-scale of GRB shocks. Hence, the 
growth of the field does not have observational conse- 
quences. Similarly, the typical correlation length of the 
magnetic field is much smaller than the source size and 
cannot be resolved. Thus, conventional light-curve obser- 
vations are unable to test the magnetic instability mecha- 
nism. However, polarization measurements might be more 
promising, as we show next. 

4.1. General Considerations 

Synchrotron radiation produced by relativistic electrons 
is known to be highly polarized, predominantly in the di- 
rection perpendicular to the local magnetic field ( Ginzburg 
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1989; Rybicki fc Lightman 1979 ). It was shown in §|| that 
the generated magnetic field is randomly oriented in the 
plane of the shock front. 

The afterglow radiation emitted by any infinitesimal sec- 
tion of the GRB blast wave is relativistically beamed to 
within an opening angle 6'b ~ 7s^^ ^ 1. Hence, an exter- 
nal observer sees a conical section of the fireball, as defined 
by this opening angle. In addition, the rapid deceleration 
of the fireball reduces its surface brightness as it expands. 
For a 



statistical properties of the ISM turbulence. If, however, 
the source is extended, then the overall pattern is obtained 
from the superposition of the incoherent patterns of its in- 
dividual parts. Thus, if the angular size of the source, 
Os, is larger than the characteristic scale of the speckles, 
namely 9s > Oq, then the intensity fiuctuations wash-out 
and the scintillation amplitude diminishes. The observa- 
tions of a late-time decline in the amplitude of i ntensity 



particular observed time, emission along the line- al. 1997 



scintillati o ns for the radio afterglo ws GRB970508 ( Frail el 
IWaxman, et al. 1998|) and GRB980329 ([Tay 



of-sig] it axis to the source center suffers from the shortest 



geometric time-delay, and hence originates at a larger ra- 
dius and is dimmer than slightly off-axis emission. The 
source therefore appears as a narrow limb-brightened ring 
(Waxman 1997c; Sari 1998; Panaitescu & Meszaros 1998; 
Granot et al. 1998a). The outer cut-off of the ring is set 
by the sharp decline in the relativistic beaming at angles 
greater than . Interestingly, the shock surface appears 
to a distant observer as almost perfectly aligned along the 
line-of-sight at the edge of the ring. This effect results 
from relativistic aberration (Rybicky & Lightman 1979, p. 
110), i.e. the Lorentz transformation of angles from the 
shock frame (in which the normal to the shock surface is 
inclined at an angle 7^"^^ relative to the line-of-sight) to the 
observer frame. Therefore, at the limb-brightened edge of 
the ring, the small-scale magnetic field is oriented tangen- 
tially on the sky. Consequently, the random magnetic field 
does not average-out but rather produces linear polariza- 
tion which is oriented radially from the center at any point 
on the ring. The resulting synchrotron radiation obtains 
a deg ree of polarization of 



Pe 



1 



Pe + 7/3 



72% 



(12) 



for the typical value of the power-law index, pe — 2.5, of 
the electron energy distribution, dNe/d'^e oc Je^', in GRB 
sources. 

The two-stream mechanism for the amplification of the 
magnetic field can be tested only if the source is resolved, 
since the net polarization of a circularly-symmetric image 
is zero^ . There are two ways for resolving a compact GRB 
source: (i) scintillations of radio afterglows due to elec- 
tron density irregularities in the ISM of the Milky Way 
galaxy (Goodman 1997); and (ii) gravitational microlens- 
ing due to an intervening star along the line-of-sight (Loeb 
& Perna 1998). Since lensing occurs only rarely, we focus 
our discussion on the first method. Observations of in- 
terstellar scintillations probe angular scales of order a few 
micro-arcseconds ( /ias), far below the VLBI resolution 
(- 300 //as). 

The interstellar scintillations arise when fiuctuations in 
the electron density randomly modulate the refractive in- 
dex of the turbulent ISM. As a result of random focusing 
and diffraction of the electromagnetic wave, a point source 
produces a spatial pattern of random bright and dim spots 
— the speckle pattern. The source brightness fiuctuates as 
the observer moves across the pattern. The characteristic 
angular correlation length of the pattern, Oq , is set by the 
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Rs ^ 10^' cm at times of ^ 1 month and ^ 2 weeks after 
these bursts, respectively. These estimates are consistent 
with the simplest fireball model predictions. 

The Weibel instability mechanism predicts that differ- 
ent segments of the ring-like source emit synchrotron ra- 
diation which is linearly polarized along the radial axis, so 
that the net polarization vanishes when averaged over the 
source. If ds ^ 6*0, the source is effectively point-like and 
hence symmetric. This regime is characterized by strong 
intensity scintillations and weak polarization fluctuations. 
In contrast, when 9s > 9o, different parts of the source 
are mapped differently, and the source is resolved. As the 
earth moves through the scintillation pattern, an observer 
will measure fluctuations in the direction and amplitude 
of the polarization, while the intensity would vary only 
weakly due to the overlap of the separate speckle patterns. 
The polarization scintillations should therefore be strong 
when the flux fiuctuations are weak. 

We consid er two types of scintillations, d iffractive and 



1998) provide an estimate for the shock radius, 
1 month and 



refractive (Goodman fc Narayan 1985; Blandford & 



Narayan 1985| ). Diffractive scintillations occur when the 
source is nearly point-like, 9s <ti 9^, relative to 



10 GHz 



-11/5 



/ias, 



(13) 



which is the diffraction a ngle for a typical scattering mea- 
sure of 10~^-^m~-^°/^ kpc (Goodman 1997). The flux mod- 
ulation amplitude in the strong scattering regime is close 
to 100%. For a Kolmogorov spectrum of ISM turbulence, 
the characteristic speckle length is 



9n ~ 2.3 I 



-lOGHz 



y/5 



/ias. 



(14) 



assuming a scattering screen distance of ^ 1 kpc and a 
typical scattering measure of lO~^-^m.~^^/^ kpc (Goodman 
1997). The time-scale for diffractive scintillations is 



idiff 



(10 GHz) 



6/5 



hr 



(15) 



if the transverse velocity of the line of sight is dominated 
by the earth with v± ~ 30 km s~^. As long as 9s <C ^Oi 
the polarization is close to zero, but when the source ap- 
proaches the speckle correlation length, 9s ~ 6*0, the po- 
larization scintillations could grow up to a large ampli- 
tude, of order a few tens of percents [cf. Eq. ([l2|)]. For 
these scintillations to be detected, the source must be 



net polarization signal might still result from an asymmetric source (e.g., due to a misaligned jet) or due to an inverse cascade of the 
magnetic field to large scales (Gruzinov & Waxman 1998). 

^"Effects H11P tn Hifi^prpntial Faraday rotation or anisotropy of the ISM turbulence are unimportant because of the smallness of the scattering 
angle, ~ /ias (Narayan 1999). 
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observed at relatively low frequencies (Goodman 1997), 
namely v <^ lOGHz for typical ISM conditions. Unfor- 
tunately, the synchrotron self-absorption often occurs at 
frequencies below 5 GHz, and so the afterglow might be 
fainter at these low frequencies, making the detection of 
polarization scintillations more difhcult. In addition, the 
source image resembles more a filled disk rather than a 
hollow ring at low frequencies (Granot et al. 1998a). The 
unpolarized radiation emitted near the center of the disk 
will thus lower the overall degree of polarization. 

As the source gets larger, 6s ^ Od, the diffractive ef- 
fect weakens, and the scintillations are dominated by the 
refractive effect, which yields only modest intensity fluc- 
tuations with an amplitude ^ 10%. The polarization fluc- 
tuations in this regime have a corresponding amplitude of 
only a few percents. The characteristic time-scale for the 
refractive modulation is 



rcf 



14 



7off 



10 ^as 



hr, 



(16) 



where 6eS is the effective size of the source (see poodmaij 
1997[]or details). 

4.2. Polarization Scintillations: Formalism 

The properties of the radiation field are fully de- 
scribe d by four scala r parameters — the Stokes param- 
eters ( iGinzburg 1989| ), which are additive for incoherent 
sources. For synchrotron radiation produced by relativis- 
tic electrons, these parameters include the intensity / and 



q being the spatial wave-number. The quantity is the 
phase structure function which yields the phase shift along 
different paths and is determined by the ISM turbulence. 
The inferred value of /? for the Galactic ISM is somewhat 
uncertain but close to the Kolm ogorov theory prediction 
f3 = 11/3 (Armstrong et al. 1995). In calculating the scin- 
tillation indexes below, we adopt the approximate value of 
/3 « 4 for which the W is Gaussian, which greatly simplifies 
the calculation. 

The Fourier transform of W is the apparent brightness 
distribution of the scattered image of a point source: 



W{e,cb)^{lo-Io)/Io 



w, 



(20) 



where ^ denotes a Fourier conjugated pair, and 9 = 
r / const and (j) are the radial and angular polar coordi- 
nates on the sky relative to the source center. 

The scattered image of an extended source is the con- 
volution of the image kernel of a point source with the 
brightness distribution at the source, Pj{d,(j)), 

w{e -0',(p- <j)') Pi{e\ (j)') 9'de' d^'. 



(21a) 

Similarly, the "images" of the other Stokes parameters are 
Qi9,^) ^Pq, U{9,dp)^W *Pu. (21b) 



Q = I cos 2ip cos 2x, 
U — I cos 2-0 sin 2%, 
V = Q. 



(17a) 
(17b) 
(17c) 



The last parameter, V , describes circular polarization 
while Q and U describe linear polarization, x is the angle 
between the polarization axis and an arbitrary fixed direc- 
tion in the sky and cos2ip = — /_l)/(/|| + /_l) is the 
difference between the radiation intensity along the two 
orthogonal axe s of polarization divided by the sum (see 
Ginzburg 19"89| ; and [Rybicki fc Lightman 1979| , p. 180). 
Both x(r) and V'(r) are determined by the source, but are 
not affected by the scintillations. The degree of polariza- 
tion is defined as 



{Q^ + U^ + V^f^ II 



(18) 



Given a power spectrum of electron density fluctuations 
in the ISM, the statistics of speckles in a scintillation pat- 
tern is usually characterized by the second moment corre- 
lation of the complex electric field of the electromagnetic 
radiation, 



W(Ax) = E{yi)E*{y. + Ax) 
cx exp[-i:)^(Ax)/2] 
cx exp [-const x (|Ax|)''"^] 



(19) 



where x and Ax are two-dimensional vectors on the plane 
normal to the line of sight, the "bar" denotes an ensemble 
average, and (3 is the power-law index of the power spec- 
trum of electron density fluctuations, |<5ne(g)p (x with 



Finally, the amplitude of the intensity fluctuations due 
to scintillations is determined by the so-called scintillation 
index: 

(J2^ X 1/2 



Si = 



(22) 



with analogous definitions for the indexes of the other 
Stokes parameters Sq and Su. We use angular brackets to 
denote integrals of the form, {W'^) = j[W{9,(j3)f 9<19 d4>. 
The normalized amplitude of the polarization scintillations 
is described by the scintillation indexes of the polarization 
signal Squ and the degree of polarization 5^, 



Squ = {Sq + S^) ^ 



Stt — Squ / Si. 



(23) 



4.2.1. Polarization Scintillations of GRB Afterglows 

To illustrate the qualitative properties of the polar- 
ization scintillations in GRB afterglows we consider a 
crude model for the source that simplifies the related inte- 
grals considerably. We approximate the circular source 
as having a uniform surface brightness over the region 
< 9 < 9s{t), on the sky. We also normalize the total 
flux to unity at all times since it enters only as a mul- 
tiplicative factor to the polarization indexes. The linear 
polarization is oriented along the radial direction, so that 
the polarization angle is equal to the polar angle % = (/) in 
equations ([TtI), and the degree of polarization is assumed 
to be constant over the source, tTs = 0.72 [cf. Eq. (|l^]. 
Much of the radiation from the ring-like image of a real 
source acquires this polarization level, although the overall 
polarization is somewhat degraded by emission from the 
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central part of the ring. Our estimates should therefore 
be regarded as an upper limit on the measurable polariza- 
tion amplitude. The brightness distribution function for 
the scattered image of a point source, W, is taken to be a 
Gaussian with a variance set by the speckle angular scale, 
W = exp[— 6'^/6'§]. The angular size of the source as a 
function of time, Os{t), was evaluated by Waxman et al. 
(1998). For a cosmological source at a redshift Zg ^ 1, it 
reads 




(24) 

where E is the total energy of the fireball and t is elapsed 
time from the detection of the explosion. The scintillation 
indexes can then be numerically calculated as functions of 
0s{t)/0o, using equations (|l^-(||). 

The temporal evolution of the scintillation indexes for a 
source with Zg = I, E = 10^^ ergs and nisM — 1 cm"'^ is 
presented in Figure ^. At early times, when the source size 
is small {9s <^ 9o), the polarization fluctuations are weak 
while the intensity fluctuations are at maximum. When 
the source size approaches the diffractive scattering angle, 
9d, the source is resolved and the observed radiation is par- 
tially polarized. At the same time, the intensity fluctua- 
tion amplitude declines due to the overlap between speck- 
les. The polarization fluctuations peak when 9s ~ at a 
value of ~ 20% x (7rs/0.72). As the source size increases 
even further, the fluctuation amplitude of both the inten- 
sity (Si) and the polarization (Squ) decrease, due to the 
overlap of scattering patterns from different regions of the 
source. However, the fluctuation level of the degree of po- 
larization (iStt) continues to increase with increasing source 
size and asymptotes at ~ tTs = 72%. Thus, the saturation 
level of is independent of the details of the scattering 
processes and provides information about the intrinsic de- 
gree of polarization at the source. 

5. CONCLUSIONS 

We have shown that the relativistic two-stream mag- 
netic instability is capable of producing strong magnetic 
fields in the internal and external shocks of GRB sources. 
The generated fields are randomly oriented in the plane of 
the collisionless shock front, and fluctuate on scales much 
smaller than the size of the emission region. The instabil- 



ity inevitably produces magnetic fields with the magnetic 
energy parameter of ~ 10^^-10^"* due to the isotropiza- 
tion of the electrons at the shock (see, e.g., the simulations 
by Kazimura et al. 1998), and could saturate at yet higher 
values of es ^ 0.1 if the protons do the same. Numerical 
simulation of electron-proton plasmas are necessary in or- 
der to examine under which conditions the protons might 
enhance the magnetic energy up to these high values. 

Galama et al. (1999) suggested a distinction between 
two classes of GRB afterglows: radio-weak GRBs like 
GRB971214 or GRB990123 where the magnetic energy pa- 
rameter might be as low as es ^ 10~^-10~^, and radio- 
loud GRBs like GRB970508 where es ^ 10~^ (Waxman 
1997a,b; Wijers & Galama 1998; Sari et al. 1998). Low- 
field afterglows are short and dim in the radio (and account 
for the majority of the afterglow population) while high- 
field afterglows are long-lived and bright in the radio. In 
our model, low-field GRBs would arise naturally due to 
the saturation of the instability at the initial kinetic en- 
ergy of the electrons. High-field afterglows might result 
from proton amplification of the magnetic energy. 

Our model for the magnetic field generation predicts 
the existence of polarization scintillations in the radio af- 
terglows of GRBs. Since the typical correlation length of 
the generated magnetic field is very small, no net polar- 
ization is expected in the absence of scintillations, unless 
the circular symmetry of the source is broken (e.g. due to 
a jet which is misaligned with the line-of-sight) or if there 
is an inverse cascade of the generated magnetic field to 
much larger scales. In the absence of such complications, 
the polarization scintillations should appear typically af- 
ter a week, when the angular size of the source becomes of 
order a micro-arcsecond, or equivalently when its physical 
size is ^ 10^^ cm. The normalized amplitude of the polar- 
ization scintillation signal at that time could be as high as 
- 10-20%. 
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APPENDIX 



ULTRA-RELATIVISTIC TREATMENT OF THE MAGNETIC INSTABILITY 
Starting with the kinetic equation 

dtf + v-d^f+{e/c)vxB-dpf = 0, 



(Al) 



for the collisionless plasma, separatin g the PDF into an unper turbed part and an infinitesimal perturbation, / = F{p) + f, 
and specifying -F'(p), one can obtain (Yoon & Davidson 1987) the following dispersion relation for the magnetic (Weibel) 
instability in the relativistic regime: 



PI 



2(l-/3i) 



(A2) 



where /3|| = p\\/^mc, l3± = p^/^mc, 7 = (1 - /3| 



(2/3^)-i ln[(l + /3^)/(l - /3^)], and py abd pj 



the components of particle momentum averaged over the PDF. Here we denote quantities parallel and perpendicular with 
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respect to the direction of the shock propagation, opposite to the convention used by Yoon fc Davidson 1987 . It is easy 
to demonstrate that the instabihty occurs for the range of k'^ given by 



0<k^ <kt 




2/31(1 



G{I3±) 



(A3) 



and only with anisotropic PDFs for which the expression in square brackets is positive. 

The mode with the largest growth rate dominates in the evolution. We therefore want to find the maximum growth 
rate, Fniax, and the corresponding wave vector of the fastest growing mode, fcmax- Upon straightforward but lengthy 
calculations, we obtain: 



k^ = 



1-/3 



2 i 2/3iG(/3i 



(1 - /3i)3/2 [l-^l 



111- 



+ (1 - 2/3i - G(/3i 



7c2(l-/32) 



-PI 



2(l-/3i) 



- G(/3±) 



(l + /3i)/3|| 

V2(l - I 1 - /3i 



1 - 2/3i - /3i 

/3i 



1/2- 



G(/3j 



(A4a) 



(A4b) 



These exact equations may be greatly simplified by assuming that the plasma is ultra-relativistic and the particle parallel 
momenta (associated with the bulk motion) are much larger than their perpendicular ones (due to their thermal motion): 
7|[ S> 7± 3> 1. Then 7 ~ 7|| = 7, and we readily obtain. 



7 



1 



2V22i 

7 



1 



V2 7i 



1 



3 7^ 
\/2 7 



(A5) 



Note that in the second equation, cjp is divided by 7^, which is much smaller than 7. 

ASYMMETRIC TWO-STREAM INSTABILITY 

Cold beam - Plasma Instability 

Here we consider the case when two interpenetrating coUisionless plasma streams have different densities and speeds in 
the center of mass frame. Instabilities which occur in such a situation are often referred to as beam - plasma instabilities. 
The lack of symmetry in the system complicates analytical, fully relativistic analysis and requires numerical simulations. 
Below we provide quantitative estimates based on extrapolation of the nonrelativistic results to the ultra-relativistic case. 



The non-relativistic case of a beam-plasma instability has been considered in different regimes (see e.g., Akhiezer et al 



1975). If the densities of the two streams are very different from each other, the center of mass frame coincides with the 



rest frame of the denser stream, which we refer to as the "bulk plasma." The lower density stream is moving with some 
velocity u relative to it and is referred to as "beam" . We denote the p arameters of the bea m by a prime. The dispersion 
relation for the magnetic instability in the case of a cold beam reads (Akhiezer et al. 1975, v.l, p. 306) 



2 12 



(Bl) 



where u is the beam velocity. We can then find the maximal growth rate and the fastest growing mode, as in Appendix 



Tmax = ^maxC^ - Wp^Wp^ (w/wthe) • 

This result suggests the following scalings with the density ratio of the beams 

Tmax OC fcmax OC {n'^/Ue)^^^ , €3 OC [n'^/rief^'^ . 

Hot beam - Plasma Instability 
When particle pitch-angle scattering at a shock is strong, the beam becomes "hot", u 



dispersion relation becomes (Akhiezer et al. 1975) 



(B2) 
(B3) 

^ uth • Then, the 



(B4) 



10 



The instability occurs when 



"the 



12 



Thus, the instabihty shuts off lov fc — )■ when 



(B5) 



(B6) 



which is satisfied when n' /n < 1. 

In this case, however, Langmuir (longitudinal, electrostatic, high-frequency) waves are efficiently generated with the 
(maximum) growth rate comparable to that of magnetic instability in the previous cases: 



Langmuir — 24/3 



1/3 



(B7) 



Random electric fields of Langmuir turbulence scatter plasma particles and provide effective collisions at the shock, so 
that the MHD approximation is applicable. A detailed analysis of this process is, however, beyond the scope of this paper. 
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Fig. 1. — Illustration of the instability. A magnetic field perturbation deflects electron motion along the a;-axis, and results in current sheets 
(j) of opposite signs in regions I and II, which in turn amplify the perturbation. The amplified field lies in the plane perpendicular to the 
original electron motion. 
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Fig. 2. — Scintillation indexes of (i) intensity, Sj, (ii) the polarization signal, Squ, and (iii) the degree of polarization, 3,^, as functions of 
time for E = 10*^ ergs, Zs = 1, niSM = 1 cm~^, and tTb = 0.72. The thin vertical line marks the time when Os = Oo- 



